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Abstract 

Wc demonstrate that // O e slcipton mixing is signific;antly more restricted than previously 
thought within the already remarkably flavor-safe i?-symmetric supersymmetric standard model. 
We calculate bounds from /x — >■ ej, /x ^ 3e and, most importantly, /x — >■ e conversion. The 
process /i — > e conversion is significantly more restrictive in i?-symmctric models since this 
process can occur through operators that do not require a chirality-flip. We delineate the 
allowed parameter space, demonstrating that maximal mixing is rarely possible with weak scale 
superpartners, while 0(0. 1) mixing is permitted within most of the space. The best approach 
to find or rule out /x O e mixing in i?-symmetric supersymmetric models is a multi-pronged 
attack looking at both /x — >■ e conversion as well as /x ^ ej. The redundancy eliminates much 
of the parameter space where one process, but not both processes, contain amplitudes that 
accidentally destructively interfere. We briefly discuss implications for searches of slepton flavor 
violation at the LHC. 



1 Introduction 



Lepton flavor violation (LFV) is predicted to occur at an unobservably small rate in the Stan- 
dard Model (SM). In low energy supersymmetric theories, new sources of lepton flavor violation 
are generic in the soft breaking sector. The experimental non-observation of — )• e processes is 
particularly restrictive, given the impressive bounds on — )• 67 from MEGA |lj and MEG [2]; on 
/i — )• e conversion from SINDRUM II [3], and to a lesser extent from — )• 3e from SINDRUM [Ij. 
Further progress is expected from the varied experiments that are ongoing as well as planned future 
experiments such as Mu2e [5J and other proposals utilizing Project X at Fermilab 

In the minimal supersymmetric standard model (MSSM), -f-)- e mixing is severely constrained 
by these bounds (e.g. 13 |8l [9l [lOl |TT] ) . The size of the mixing can be characterized by the quantity 
S^Y = (5m^y/m^ where (5m^y is the off-diagonal (12)-entry appearing in the sfermion mass matrix 
connecting the X-handed slepton to the y-handed slepton, and is the average slepton mass. 
Ref. [m found 10-^ while < 6 x 10"^ over a scan of the mSUGRA parameter space. 

Similarly strong bounds on 6^j^j^ can also be found, though cancellations between diagrams in the 
amplitude can in some cases allow for much larger mixing [9l [TU[ [TT] . 

Recently, a new approach to weak scale supersymmetry that incorporates an extended R- 
symmetry [T2], suggests large flavor violation in the supersymmetry breaking parameters may 
be present without exceeding the flavor-violating bounds. This is possible for several reasons: R- 
symmetric supersymmetry has no flavor-violating LR mixing, solving the worst of the problem 
trivially. i?-symmetric supersymmetry has Dirac gauginos, and no Majorana masses, removing 
all dimension-5 flavor-violating operators. Finally ii-symmetric supersymmetry also has no flavor- 
conserving LR mixing, and so there are no "large tan (3 enhanced" effects. These beneflts were 
found to virtually eliminate constraints on the slepton flavor mixing |12j . 



In this paper we reconsider the constraints on slepton mixing, speciflcally, fx -H- e mixing. 
Unlike the MSSM, the most important constraint is not necessarily fi — )• ej. This is easily seen 
by inspection of the /^-symmetric flavor-violating operators: ^ — )• 67 requires a chirality-flip via a 
muon Yukawa coupling, whereas /U — )• e conversion has no such requirement. We flnd that /i — )• e 
conversion rules out maximal mixing throughout the right-handed slepton mixing parameter space 
for sub-TeV superpartner masses. This is complementary to — >• ej, where we find cancellations 
between the bino and Higgsino diagrams, analogous to what was found before in the MSSM [9l 
[TOl lllj. For left-handed slepton mixing, we find possible cancellations in the amplitudes for — )• e 
conversion, and instead ;U — t- 67 provides generally the strongest constraint. We also calculated 
^ — 7- 3e and find it provides the weakest constraint on both left-handed and right-handed slepton 
mixing throughout the parameter space we consider. 

This paper is organized in as follows: We review the relevant characteristics of a model with an 
extended ii-symmetry, and the super GIM mechanism in Sec. [2} In Sec. [3j we begin the discussion 



of experimental constrains on the parameters from /i — )• 67, in Sec. |3.1[ /x — )• e conversion in Sec. 3.2 
and finally, /i — t- 3e in Sec. |3.3[ In Sec. [4] we briefiy discuss implications for slepton flavor violation 
to be observed at LHC. Finally, in Sec. [5] we conclude with a discussion of our results. 
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2 A Simplified R-symmetric Model 



We are interested in analyzing LFV in the minimal i?-symmetric standard model (MRSSM). The 
gaugino structure of the MRSSM has been studied in detail in Ref. [13], where the mixings and 
couplings of the four Dirac neutralinos and four Dirac charginos are given. Weak scale supersym- 
metry with Dirac gauginos is a possibility that was contemplated some time ago [HI [151 US] and 
more recently [I71[I3[I3[2ni[2Il[23[23l[I2l[Il[M[23[2ni[Z3[2Hl[Ml[^ A fully general 

analysis of LFV in the MRSSM would be a substantial undertaking. Fortunately, there are sev- 
eral simplifications we can employ to gain a fairly general understanding of the allowed parameter 
space of LFV in the MRSSM. One important restriction is that the Dirac wino cannot be light 
in the MRSSM, due to the structure of the wino supersoft operator |17J. Essentially there is an 
unavoidable contribution to the vev of the 5C/(2)2,-triplet scalar that causes a contribution to the p 
parameter that is too large unless the wino is above about a TeV. Secondly, since there is no cou- 
pling between up-type Higgs and leptons, the contribution from the up-type Higgsino eigenstates 
is suppressed by the small mixing between bino or Hd and H^, and so can be ignored. 

Itemizing the simplifications, we take: 



1. The wino mass, M2, is taken to be sufficiently large so as to give negligible contribution to 
flavor violating interactions. This simplification means that the p-parameter is automatically 
safe throughout the parameter space we consider. 

2. The up-type Higgsino mass /i^, is also taken to be large for convenience. Since the up-type 
Higgsinos play no role whatsoever in charged lepton flavor- violation (given also point 1), this 
is done simply to keep the gaugino sector to a 2 x 2 structure and thus easily understood. 
(We will, however, consider effects of a light up-type Higgsino on flavor-violating signals at 
LHC in Sec. |4|) 

3. We consider left-handed and right-handed slepton mixing separately. This is standard prac- 
tice when considering flavor- violation in the MSSM (e.g., [H]). We will see that there are 
qualitative differences between the allowed parameter space of left-handed and right-handed 
slepton mixing. 

4. We assume the slepton mixing is purely in the 2x2 flavor space of e, /i. Enlarging this mixing 
to the full 3x3 mixing does not qualitatively change any of our results, and instead simply 
dilutes the effect of the mixing, while adding more mixing angles and thus more parameters 
to the model. Since the focus of this paper is to explore e mixing, no further discussion 
of the 3x3 case will be given. 

5. For our numerical results, we take m^-^ = 1.5m^-^. This seems a far more drastic assumption 
than it actually is. Our motivation is to consider slepton flavor violation when there is 
essentially no degeneracy among the sleptons, and so we took the slepton mass ratio to be 
"order one" but not near one. Taking the ratio much larger than one does not appreciably 
increase the flavor violation, while taking it smaller causes the super-GIM mechanism to 
suppress the flavor-violating signal. Our compromise is the above number. 
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In Appendix [XJ we provide more details on the gaugino structure and flavor- violating interac- 
tions as directly relevant to this paper. With the above assumptions, there is only one light Dirac 
chargino (which is H^-like) and two light Dirac neutralinos (which are mixtures of and B). 

A few more comments on the slepton mass eigenstate hierarchy are in order. MSSM analyses 
of slepton flavor violation have, by necessity of LFV constraints, concentrated on the case where 

the mass difference between the different states is small, = m? — m? <C m? . In this limit, it 

h h '1,2 

is straightforward to show that the contribution to LFV can be expanded in powers of A^, taking 
the form 

sm2e^(-^ + ...) , (1) 

where MsusY is typically the largest mass sparticle in the diagram that dominates the process. 
There is no A-independent contribution within the parentheses due to the super GIM mechanism 
(see the next section). Since sin29i = 2mg^/A^, one factor of A^ cancels, giving proportionality 
to the 5 parameter mentioned in the introduction and used in many other papers on LFV in the 
MSSM (at least up to a possible further suppression of \m^^mjij / rnj^ if mN,c ^ "^[j 

In this paper, A^ is not small, and so using the "(5 parameter" is simply not appropriate. Instead, 
it is easy to see that in the opposite limit, A^ ^ m? , the contributions to LFV are proportional 

to simply sin2^;/m? . Hence, the relevant parameters we show in most of our numerical results 
are bounds on sin 26*; as a function of the slepton, gaugino, and Higgsino masses. Reducing the 
splitting can be roughly approximated by relaxing the constraint on sin 2^; by ratios of A^j^/A^^^. 



2.1 The super GIM mechanism 

The "super-GIM mechanism" - the GIM mechanism applied to flavor in the superpartner sector - 
is important in understanding the phenomena of flavor violation. As is well known, the super-GIM 
mechanism arises as a consequence of the unitarity of the slepton mixing matrices that diagonalize 
the mass matrix; ul^Ukj = where the sum over repeated indices is performed. This combination 
of mixing matrix elements always appears as a prefactor in the calculation of amplitudes of flavor 
violating processes. Specifically for our case of slepton flavor violation, we have ulf^Uk/j. = 0, 
corresponding to an incoming muon, and an outgoing electron, with internal sleptons labeled by 
k. The sum over k corresponds to summing over all mass eigenstate sleptons Ik in the loop. There 
are two immediate consequences of the super-GIM mechanism. 

First, terms that do not depend on the slepton masses do not contribute. Let /(m^) be some 
function that depends on the mass of the sleptons and a be some quantity that does not depend 
on nik, then 

^ UlUk^[a + f{mk)] = UlUk,f{mk). (2) 

k k 

The form of Eq. ^ appears when a logarithmic divergent loop integral is dimensionally regularized, 
and one finds the 1/e term appearing as a constant term a in the above equation. This leads to 
an important result: the would-be logarithmic UV divergence in flavor-conserving processes is, in 
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fact, UV finite in flavor-violating processes. In this paper, unless otherwise stated, we will omit the 
terms in our expressions that are canceled by the super- GIM mechanism. 

The other well known consequence is that, when all the sleptons are degenerate, there is no 
flavor violation. This can be seen again in Eq. ^ with rrik = m, the sum over all slepton flavors 
in a flavor-violating process vanishes. 



3 Experimental constraints 

There are three — t- e conversion processes with experimental bounds: /i — t- ej, // — )• e conversion, 
and /U — 7- 3e. In this section we present our calculations of the rates of these processes and present 
results in terms of a series of contour plots showing the allowed parameter space. 

The rate for // — t- 67 was estimated in Ref. |12| in the slepton flavor-violating mass-insertion 
approximation with a pure bino and wino and a specific gaugino hierarchy. In this paper we have 
neglected the wino, due to the p parameter constraint, and instead included the down- type Higgsino 
H^. Since we have considered large mixing angles, up to and including maximal mixing, we have 
diagonalized the slepton masses explicitly and done our loop calculations involving the slepton mass 
eigenstates. 

As stated in our simplifications, we have not included contributions from the wino or up-type 
Higgsino. We focus on the case where the sleptons and the lighter neutralinos are in the sub-TeV 
range where wino contributions can be reasonably ignored. The up- type Higgsino does not couple 
to leptons, and we take the light quark Yukawa couplings to vanish. Thus, the up-type Higgsino 
does not give a significant contribution to any of /i — t- 67, /i — t- e conversion in nuclei and — >• 3e. 

With these simplifications, the amplitudes of LFV processes are sensitive to just two neutralinos, 
mixtures of B and inside the loops. We can also neglect the contributions due to charginos 
because the only light chargino is H^-like. Hence, all types of diagrams we consider involving a 
chargino are suppressed not only by one power of muon Yukawa, but also one power of either the 
electron Yukawa, or the tiny wino content of the light chargino at the lepton-chargino-sneutrino 
vertex. This also means that sneutrino mixing does not contribute to LFV processes, and thus the 
difference in the amplitudes between left-handed and right-handed slepton mixing is due solely to 
the hypercharges and masses of the left-handed and right-handed charged leptons. 

3.1 /i — )■ 67 

The neutrinoless muon decay /i — >■ 67 occurs through the effective magnetic dipole moment operator, 
ea^iyF^'^ fi, and requires a chirality flip of fermions. There are no tree level operators that lead to this 
decay, and the lowest order is at one loop. In the MRSSM, there are only two types of contributions 
to the /i — )• 67 amplitude: one where the chirality flip occurs on the external muon line, and the 
other where the flip occurs as a result of a muon-smuon-Higgsino vertex proportional to the muon 
Yukawa coupling. The diagrams are shown in Fig. [T| 
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(a) Chirality flip on the external muon line. (b) Chirality flip at the Yukawa vertex. 




Figure 1: Feynman diagrams for — > 67 corresponding to the amplitudes (a) Ajuni and (b) Ajiin2 
mediated by right-handed slepton flavor mixing. The diagrams for left-handed slepton flavor mixing 
are obtained by swapping L R. 



We calculated the amplitudes in the mass eigenstate basis of the sleptons and neutralinos, and 
as a check we derived the results obtained in Ref. [8] (replacing their Jl-f mixing with e-/i mixing). 
The effective Lagrangian is 



m 



eff 



(3) 



We rewrite the amplitudes ^^^^ip ^^'^ ^^dm^ 



7dip' 



-^7dip 



7dip 



i=l 
2 



(4) 
(5) 



i=l 



where the sum is over the i-th neutralinos. The subscripts 1 and 2 denote the locations of the 
chirality flip on the muon line and at the muon-slepton-gaugino vertex, respectively. As we shall 
see below, for right-handed sleptons there can be an accidental cancellation between amplitudes 
involving these diagrams. 



The /i — 7- 67 branching ratio is given by 



BR{n 67) 



|4i |2_L|zl« 
I 7aipl ^ l^-vr 



7dip I 



(6) 



with the amplitudes involving a neutralino Ni and sleptons li and I2 with the sleptons mass-ordered 

ided sleptons are 

fnl{xii) fnl{x2i) 



as m^^ < THj^. The amplitudes involving right-handed sleptons are 



A 



Rinl 



2 

rrir 



rrir 



Yy^ruNi 



A 



Rin2 



2(167r2)Mzsin6l^ cos/3 



OmmOLiB cos 6 1 sin 9i 



fn2{xii) fn2{x2i) 



rrir 



(7) 

(8) 
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where Apuni is the amphtude that involves an external chirality flip of the muon and Aiiin2 involves 
a flip at the Higgsino vertex. Here O^-f^o and O^-jj are the Higgsino and bino content of Ni^ re- 
spectively (i.e., the corresponding elements in the orthogonal matrices that diagonalize the gaugino 
mass matrix squared), and = Y^" = +1. To lowest non- vanishing order in M^, the neutralino 
mixings are (dropping the subscripts L and R from now on): 

r, ( i^/i \ n ( ^ i\/r \ cos l3 sin B^Mz fid 

OiB(^d < Ml) = Oi^o(/id > Ml) = —2 2 ' (9) 

^ / , ^ N ^ / , ^ X COS 13 sin OrnMzMi , , 

02^(;Ud»Mi) = 02^„(/x,«Mi) = \ 2 \ (10) 

and 0-(r, ttO) = 1 in the appropriate limits. The functions /„j(xi), with = /m? , with 
j = 1, 2, come from integrating over the loops in the diagrams: 

fniix) = —^-—j{l-6x + 3x^ + 2x^-6xHnx), (11) 
2(1 — xp 

fn2{x) = —^—.{l-x^ + 2x\nx). (12) 
(1 — x)^ 

Finally, the amplitudes for the left-handed sleptons can be obtained from the right-handed slepton 
results by doing the replacements 

A^^^^Ai;^^ upon {Ylm}j^{Ylrrqj. (13) 

Inserting the results in Eqs. ([9])-(10) into ([7])-([8]), we see that to lowest vanishing order in Mz, 
BR{^ — )■ ej) is independent of tan /3. We can also see explicitly that when the two slepton masses 
are degenerate, the branching ratio vanishes, as expected from the super GIM mechanism. 

As an aside, it is also straightforward to see what happens to the results when the mass hierarchy 
between the slepton and the neutralino are inverted. The loop functions satisfy the identities, 

/nl(^) + /ni0) = \, (14) 
Xfn2{x)- Ul^] = 0. (15) 



We are now in a position to discuss the amplitudes in various limits. In the bino-like limit 

Ml <^ fid, one sees that Ajuni dominates, as Ajiin2 is of order Mi/ ^d- 

When A'^i becomes i?^-like, there is a cancellation between the amplitudes involving a chirality 
flip on the external muon line, and the one with the flip occurring at the muon Yukawa vertex. The 
dominant diagram in the S-like case, ^ijini, is now suppressed by fj.'^/Mf, the same suppression 
factor appears A/ji„2- So the dominant amplitudes come from the diagrams involving a i?-like 
neutralino exchange. Note that ^_R2n2 has an opposite sign compared to ^H2ni and the total 
amplitude can vanish for some choice of parameters. 

In Figs. [2]^a)-[2]^d), we show the allowed regions in MRSSM parameter space with right-handed 
slepton mixing that satisfy the bound BR{ix — ;> 67) < 1.2 x 10^^^ [H [2]- 
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The situation is drastically different in the case of left-handed slepton mixing. The hyper- 
charge of the left-handed leptons (Y^ = —1/2), has an opposite sign to the right-handed lepton 
hypercharge, and so the amplitudes interfere constructively, instead of destructively as in the case 
of right-handed slepton mixing. This leads to a more severe bound on the allowable regions in 
parameter space for left-handed slepton mixing. This is shown in Figs. [3]^a)-[3]j^d). 



3.2 — )■ e conversion in a nucleus 



The conversion of a muon into an electron can give a qualitatively distinct bound on /x -f-)- e slepton 
mixing because there are several types of operators beyond those that contribute to /i — )• 67. We 
discuss the operators for — )• e conversion, one-by-one, in this section. 

The /i — )■ e conversion amplitude is dominated by coherent processes, and so we only took the 
quark vector currents into account. The operators that contribute to the incoherent terms, q^y^q, 
q^'^j^q, and qa^'^q have been neglected. This leaves us with the scalar and vector current, qq and 
q^'^q, respectively |7]. 

The only diagram that can contribute to a scalar quark current is the box diagram. Without 
left-right mixing of sleptons in the MRSSM, the dominant term, with bino couplings at each vertex, 
contains no chirality flip of the quarks, and is therefore a vector current. We also take the quark 
current to be non-relativistic to simplify the calculation involving the magnetic dipole term. Thus, 
the amplitude for — )• e conversion is well approximated, for our purposes, by only taking quark 
vector currents into account. 

The diagrams we consider are the photon penguin, the Z penguin, and the box diagram shown 
in Figs. 4|5[6|7 We only take the dominant terms of the box and the Z penguin amplitude into 



account: that is, the terms involving the bino coupling at each vertex which does not contain any 
chirality flips of the external fermions. The effective Lagrangian at the parton level can be written 
as [7] 



q=u,d 



riA'^PL + A^Pr) + 'Ih^ta^^-kM'^A^PL + A^^,^Pn) 



+ e'Yl ^^'i^-^z + Ai'jPL + iA§ + All)Pn]f^qj,q, (16) 



q=u,a 



where Qq is the quark electric charge, fc^ ~ """^^ is the momentum transfer, A^'^ and A^'^_^ ^ 

correspond to the 7-penguin and Z-penguin, respectively, and A'^^^^^ corresponds to the box 
diagram. 

The most severe upper bound to date is on the conversion rate ratio with a gold nucleus 
BR{^jL e)Au = r(^- Au ^ e- Au)/r(^- Au)capturc < 7 x lO^^^ f^om SINDRUM II [3]. Because 
of the large number of protons in the gold nucleus, the distortion to the muon wave function from a 
plane wave must be taken into account when evaluating the overlap between the muon and nucleus 
wavefunctions. This has been done in Ref. |33], and we will use their overlap integrals, with the 
neutron density determined from pionic atom experiments (method 2 in [33j). Other nuclei could 
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Ml 

^ (a) /id = 100 GeV (b) /id = 200 GeV 




100 150 200 250 300 350 400 100 150 200 250 300 350 400 

mj^ [GeV] m,-^ [GeV] 



Figure 2: Regions in parameter space (shaded) that satisfy the /U — >• 67 bound for right-handed 
slepton mixing. The mass of the heavier slepton is set to l.bm^^. From Ught to dark, the shaded 
areas denote mixing with sm29-^ = 0.1, 0.5 and 1, respectively. The funnel regions in the plots with 
/i^ = 100, 200 GeV is caused by the cancellation between the amplitudes involving the bino-like 
and the if9-like neutralinos. 
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Ml 

^ -^Ti (a) /id = 100 GeV (b) fia = 200 GeV 




100 150 200 250 300 350 400 100 150 200 250 300 350 400 

m^^ [GeV] m^^ [GeV] 

Figure 3: Same as Fig. [2] but for left-handed slepton mixing. We have restricted Mi < 500 GeV 
since contributions from wino-hke charginos not been included (see Sec. ^ for a discussion). 
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/ / 

Figure 4: Schematic diagram illustrating the set of penguin contributions to /i — t- e conversion 
(for f = q) as well as ^ — )• 3e (for / = e). The blob in the figure arises from both charge radius 
subdiagrams shown in Fig.[5| as well as Z penguin subdiagrams, the dominant ones shown in Fig. |6j 



also be of interest, particularly as a way to distinguish different models [3lj. The conversion rate is 
r^^e = 4my\j^X + + + + {L^R), (17) 



where, 



*'7dip - --^^^dipD, (18) 

<^ = A^V^P\ (19) 

-(2^1^^, + At^^W^P^ - {A^^^ + 2At^jv'~^\ (20) 

[{2Zu + Zd)V^P^ + {2Zd + ^„)V^(")]Af , (21) 



7 7 



where = (Z,^ + ^qi)/2, with Zq^^^^^ = - Qsin^ 0^, = 1/2, = -1/2 for up and down 
type quarks, and = 0. The first term in Eq. (17), proportional to l^^jp + =f^^ + s^b^x + ■^z'l'^^ 



corresponds to slepton mixing in the right-handed sector, while the second term proportional to 
l-'^^ip + ■^■y' + -^box "I" -^z corresponds to to slepton mixing in the left-handed sector. The 
coefficients D and y^P'") are to the overlap integrals of the muon and the nucleus for the leptonic 
dipole and vector (proton, neutron) operators. We used, for a gold nucleus, D = 0.167, V^^^ = 
0.0859, = 0.108 from Ref. [33] • 

Now we will discuss each diagram below. We will present the results for both left- and right- 
handed slepton mixing. But, for simplicity, we will only discuss the case of right-handed slepton 
mixing explicitly. The amplitudes corresponding to left-handed slepton mixing can be obtained 
from the right-handed ones by replacing the appropriate hypercharges and slepton masses. Note 
that for the Z-penguin, there is also an additional minus sign after the replacement of hypercharges 
and slepton masses. 



3.2.1 Charge radius 

The charge radius amplit 
The dominant term is the one involving the U-like neutralino in the loop, with B coupling at each 



L Ft 

The charge radius amplitude A-^' comes from the 7-penguin, without a chirality flip of the leptons. 
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(a) (b) (c) 




Figure 5: Contributions to the effective vertex from the charge radius operator. Graph (c) is 
suppressed by a factor of rrig/m^ compared to (b), and can be ignored in the hmit of vanishing 
electron mass. Also in this limit, graph (b) exactly cancels graph (a) for vanishing photon momen- 
tum, satisfying the Ward identity. Only right-handed slepton flavor mixing diagrams are shown, 
while left-handed slepton flavor mixing diagrams are obtained by swapping L o i?. 



vertex connecting a lepton. The other terms are suppressed either by the muon Yukawa or by two 
powers of the small bino content in the H^-like neutralino. The contributions to the effective vertex 
of the charge radius is shown in Fig. [5] Summing over these contributions giv^ 

= -^^H^ ) - K ^ -r.)' (22) 

with 



5767r2 mf \m] 



Ux) = — ^(2 -9a; + 18x2- llx^ + 6a;^ Inx). (23) 
1 — X* 



3.2.2 Magnetic dipole 



The magnetic dipole amplitude ^^^jp is the one that appears in /i — ^ 67, which was discussed in 
detail in the last section. For right-handed slepton mixing, the amplitude of the dipole term is 

L Ft 

smaller than the charge radius term, A-y' , due to the destructive interference between amplitudes 
involving chirality flips at different locations in the diagram. The situation reverses in the case of 
left-handed slepton mixing, where both terms contributes and the magnitude becomes larger than 
the charge radius term. 



3.2.3 Z-penguin 



The Z-penguin contribution contains diagrams in Fig. [Sj with the photon replaced by the Z boson. 
The contribution coming from this set of diagrams is suppressed by 0{m'j^/M'^) compared to the 
charge radius so is negligible. Then, the dominant term is the one involving a Higgsino-Higgsino-Z 
vertex, shown in Fig. [6j 

'^We have checked that, even when = Afi, the value given by this expression differs to the exact one by < 
1%. So this expression is valid over all ranges of Mi and ^d- The discrepancy comes from the small mass splitting 
of the neutralinos when the gaugino and Higgsino masses are degenerate. We have used the exact expression in our 
numerical analysis. 
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(a) 




(b) 




Figure 6: Contributions to the effective vertex from the Z penguin. Diagram (a) gives the term 
proportional to fz in which the Z boson couples to the i?-partner of the down type Higgsino, V'f/O; 

d 

and (b) gives the term proportional to g^, with Z coupling to H^. Only right-handed slepton 
flavor mixing diagrams are shown, while left-handed slepton flavor mixing diagrams are obtained 
by swapping L R. 

We find that the Z-penguin is sub-dominant in a large region of the parameter space. The 
Z-penguin is the only amplitude that is sensitive to tan/3, and in the limit Mz <C M^, it scales as 
cos^ (3. The Z-penguin amplitude is 



sin 207 



aR . 

■g oiii 
^ ~ 647r2 Ml sin2 0^ cos^ 6^ "^'^ ' 



E 



where 

^ij = OLiiOljl 



OLi20Lj2fz 



2 ' 2 

mf mf 

h h 



'iORi20Rj2gz 



(Ml Ml 

\ 2 ' 2 

n h 



The functions fz{xi,Xj) and gz{xi,Xj) ar^ 



fz{xi,Xj) 
9z{xi,Xj) 



In Xi + 

yfXiX'j 



x: 



\n.Xi x'jlnxj 



1 Xj 1 X 'j 



Xj^ In Xj^ Xj In Xj 



1 - Xi 



1 — Xi 



(24) 



m? — ^ mr 



(25) 



(26) 
(27) 



Note that the Z-penguin effective vertex does not explicitly depend on l/Mgygy in the case 
of all other amplitudes. This corresponds to an operator of dimension-4. This is perfectly fine, 
because the weak symmetry is broken, so the weak current is not conserved. However, it is required 
that in the limit of unbroken electroweak symmetry, this effective vertex vanishes. This is easy to 
check in the limit Mz — )• 0. In this limit, the neutralinos we consider do not mix [c.f., Eq. (35)]. 



But the amplitude for the Z-penguin contain at least two powers of the neutralino mixing matrix 
elements, regardless of whether it is bino-like or Higgsino-like. Therefore this operator vanishes in 
the limit Mz — >• 0, when the electroweak symmetry is unbroken. 

^Note that the function fz appears to contain a log term that is asymmetric in the two neutralino lines in the 
loop, not as one would expect. But remember that this log term is subtracted by one containing the heavier slepton 
mass, and the final result is symmetric in the neutralinos and anti-symmetric in the sleptons, as expected. 
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Figure 7: The box Feynman diagram for — )• e conversion. Due to the conservation of i?-charges, 
the chiraUty of the squarks must be the ones shown in the diagrams. Only right-handed slepton 
flavor mixing diagrams are shown, while left-handed slepton flavor mixing diagrams are obtained 
by swapping L o i? everywhere. 



For left-handed sleptons, the Z amplitude can be obtained by replacing the appropriate hyper- 
charges and slepton masses, as well as an additional factor of (—1). This sign change arises from 
the NNZ coupling, in contrast to N'^N'^Z in the case of right-handed sleptons. 



3.2.4 Box diagram 



For the box diagram, the dominant term is the one containing bino couplings at all four vertices, 



647r2 



2^2 



where 

j4{xi,Xj,y) 



X? In Xi 



h h h 



x'j In Xj 



(28) 



+ - 



y2 Iny 



{l-Xi){xi-Xj){xi-y) {l-Xj){xi-Xj){xj-y) {1 - y){xi - y){xj - y) 



-■ (29) 



We can compare the box amplitude with A^' by approximating V^^'^ ~ V^'^\ giving 



^box 



e2 /^(x) 



3{Yi)' + 3{Y}ir] ^ 19 



u\2i 



J 4{x,x,y) 



(30) 



where x = Mf / mj and y = m| /m? . The right hand side is plotted in Fig. jsj We can see that the 
box can give a large contribution the total amplitude when the squarks are not far heavier than 
the sleptons. 



3.2.5 Numerical Results 



We took tan (3 = 3 for our analysis. The amplitudes contributing to /u — )• e conversion in gold are 
shown in Fig.|9]for right-handed slepton mixing, and in Fig. 10 for left-handed slepton mixing. The 
slepton mixing angles are taken to be maximal. For comparison, we also drew the line where the 



14 



log 




Figure 8: A plot of the right hand side of Eq. (30), 19j4(x, where x = Mf/rn^ and 

y = m^/rriy . The contours are y = 1, 10, 25 from top to bottom. The box amphtude is larger than 
the electromagnetic term when the contour is above the x-axis. 



experimental bound on the amplitude would be, as if only one amplitude were contributing to the 
conversion rate. 

For right-handed sleptons, either the charge radius or the box diagram dominate over other 
contributions. Each of these amplitudes exceeds the bound alone and they interfere constructively 
with each other. Therefore, maximal right-handed slepton mixing is excluded throughout the 
parameter space we explore. The magnetic dipole destructively interferes with the box and the 
charge radius diagrams, at small slepton masses before the magnetic dipole vanishes. However, this 
cancellation is insufficient to bring the amplitudes below the bound. 

In the left-handed slepton mixing case, the box diagram is suppressed by the left-handed quark 
hypercharge, and is much smaller. Also, the two largest amplitudes, the charge radius and the 
magnetic dipole, destructively interfere with each other, resulting in the funnel region shown in 



Fig. 12 



For both right-handed and left-handed slepton mixing, the Z-penguin is sub-dominant. More- 
over, for larger values of tan (3, the Z-penguin will be even more suppressed, since it is directly 
proportional to cos^ /3 = 1/(1 -|- tan^ /3) to lowest order, in the limit Mz ^ Afjv. We show the 



exclusion plots for fi ^ e conversion in Figs. 11 and 12 



3.3 fi^3e 

Finally, we investigate the decay fi~ — )• e^e+e^. The diagrams that contribute to this decay are 
similar to the process /i — >• e in a nucleus. While the amplitudes for this decay are not enhanced 
by nuclear factors as in the case of ^ — t- e conversion, there is a log enhancement proportional to 
log /me, arising from an infrared divergence cutoff by the electron mass. 

All of the diagrams in /x — )• 3e can be obtained from the /i — )• e conversion diagrams by replacing 
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Figure 9: The magnitudes of various amplitudes at maximal mixing of right-handed sleptons with 



degenerate squark masses of 1 TeV (i.e., the terms in Eq. (17) before taking the square). The 
contours are, .s/^ (blue), =2^^^^ (green), |=«2/^^jp| (red), and (brown). The dashed line corre- 

sponds to the bound on /i — )• e conversion as if only one amplitude were contributing. One can 
see that there are regimes where only the box and the charge radius amplitudes contribute signifi- 
cantly [subfigures (a) and (b), especially in the high Mi regions in these figures], and where all four 
amplitudes contribute significantly [subfigure (c)]. In subfigure (d), the magnetic dipole amplitude 
reaches zero near m-^_^ ~ 330 GeV. This coincides with the "funnel" region in the parameter space 
plot for // — )• 67, Fig. [2| 
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Figure 10: Same as Fig. 9 but with left-handed slepton mixing instead. The contours are, 
(blue), —ss^box (green), —My^^^p (red), and (brown). The magnetic dipole and the charge radius 
amplitudes interfere destructively with each other, opening up a large region in the parameter space 
that satisfies ^ 



e conversion. This forms the funnel regions in Fig. 12 



17 




100 150 200 250 300 350 400 100 150 200 250 300 350 400 

m^^ [GeV] m^^ [GeV] 



Figure 11: Allowable regions for — )• e conversion in a gold nucleus with right-handed slepton 
mixing. From light to dark, the shaded areas denote mixing with sin 20^- = 0.1,0.5 respectively. 
The squark masses are set to be degenerate at 1 TeV. Note that this completely rules out maximal 
mixing for right-handed sleptons in the sub- TeV range. 
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the quark line by an electron line with outgoing and e~ . All diagrams except the box are the 
same and will not be discussed here. For the box, conservation of i?-charges enforces both sleptons 
in the loop be of the same "chirality" . The box amplitude for /i — t- 3e for right-handed sleptons is 

where Ui = cos^ 6j and U2 = sin^ 9j. The factor (— 1)*"*"^ comes from the super-GIM mechanism. 
The rate for the decay /x — )• 3e is 



+ g l-DboxJ + -"box ~ 2^7dip-Obox + RR + RL 



2 4 2 8 4 

+ g-^box-^-Rfl + 3^7 -^-RiJ + 3^7 -^-RL - 3^7dip-^-Ri? - -^^^dipFRL]-, (32) 

where Frq, = A^Z^^, with a = L,R. The quantity Za is part of the electron-Z coupling; Zl = 
— 1/2 + sin^ 9^, and Zji = sin^ 0^. The branching ratio of this process is obtained by dividing the 
rate by the muon decay rate. Note that the term proportional to (^^^ip)^ enhanced by the log 
term, which is divergent in the limit of massless electrons. Our result for this divergent term agrees 
with 1351. 



In Figs. 13|14 we show the bounds on the MRSSM parameter space arising from satisfying the 



existing experimental bound BR{fi — )- 3e) < 1.0 x 10~^^ from SINDRUM [3j. The bounds on the 
MRSSM parameter space from // — )• 3e are weaker than the combined bounds from // — )■ 67 and 
// — )■ e conversion. 



4 Implications for Flavor Violation Signals at LHC 

One of the most interesting implications of the MRSSM is that flavor mixing could be at or near 
maximal throughout virtually the entire slepton and squark sector (save only perhaps for d-s 
mixing [36] )• For sleptons, this opens up the possibility of observing large ^-e mixing at colliders. 
Slepton mixing at colliders has been extensively studied p7l[551[5^linimili^li5llMII^IiBllT7l 
llHllinilSOliniESlESllMlESlES], though analyses have generally been relegated to MSSM scenarios 
where the splitting between the e,// eigenstates is very small, so as to satisfy the stringent LFV 
constraints. One of the most sensitive techniques to search for fi-e mixing is through the decay of 
a heavier neutralino to a lighter one through an on-shell slepton. This decay can arise at a large 
rate at the LHC starting with squark and/or gluino production, where the squark decays to the 
heavier neutralino and so on, such as 

q^qNi] Ni^ e^/^^[^ ; ^ ^^/e^TV,- . (33) 

The distinctive kinematic features in this cascade of 2-body decays can be utilized to extract the 
mass of the slepton through a kinematic edge (e.g. [56l [571 ISHl [SHI [6Ql l6T]). 
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Figure 13: Regions of the parameter space that satisfy the — >• 3e bound at different mixing angles 
of right-handed sleptons. The values of sin 20,- are, from hght to dark, 0.1, 0.5, 1. 
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Table 1: Mass spectrum 
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Table 2: Leading order production cross sections for squarks and gluinos at the LHC with ^/s = 
14 TeV in the MRSSM. 

In light of the bounds on the MRSSM parameter space that we have found from LFV processes, 
it is interesting to consider whether large mixing could still be seen at the LHC. A detailed collider 
study is beyond the scope of this paper, nevertheless we can use our results to uncover characteristic 
regions of parameter space where sin20; ~ 1 simultaneous with several- hundred GeV sparticles, 
and thus, where large ^ O e mixing remains within reach of the LHC. 

Closely examining Figs. [3)^d)|12[d))14[d), we discover one (small) region in the MRSSM pa- 
rameter space where the left-handed slepton mixing angle can be maximal, sin 20; = 1. For this 
region, and given first and second generation squark masses to be 1 TeV (consistent with what was 
assumed for the — )• e conversion numerical results), we compute the leading order production 
cross sections and decay rates. We take the wino mass and the right-handed slepton masses to be 
2 TeV for simplicity. The other gaugino masses in this region are Mi = 500 GeV, = 400 GeV, 
flu = 100 GeV. The mass spectrum is shown in Table [T] 

Using MADGRAPH [62], we calculated the leading order squark and gluino production cross 
sections at LHC with ^/s = 14 TeV center of mass energy for several values of the Dirac gluino 
mass for those production modes allowed by /^-symmetry in Table [2j One important observation 
made in Ref. [32j is that, for gluinos less than about 2 TeV, associated gluino-squark production 
gives the largest production rate of squarks. 

The decay rates of the squarks, neutralinos, and charginos, computed using BRIDGE [63], 
can also be computed as a function of the mixing angle 9i , shown in Table |3j For the particular 
point we considered, the first two generations of squarks decay overwhelmingly into the bino- 
like neutralino, N^. The subsequent cascade decays into opposite flavor leptons have the rates 
BR{N3 eiiNi) = O.Usin^ 26*;, BR{N3 (ee//i^)iVi) = 0.27(sin'^ 6*/ + sin^ 6*/). If the gluino mass 
is 1 TeV, for example, then the gq production leads to a total cross section of about 1 pb. With 
maximal slepton mixing, the cross section for opposite sign e// events is expected to be of order 
100 fb. Extracting this signal from background, particularly given the potentially problematic 
technique of flavor-subtraction, remains challenging. (See Ref. |[32j for a discussion of signal plus 
background analysis of a non- minimal i?-symmetric model.) 

Just as in the MSSM, one can search for the kinematic endpoint in the invariant mass distri- 
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Table 3: Decay branching ratios of the particles involved in the cascade decay N-^ — )■ Z ^ I I'^Ni 
given the MRSSM parameters given in Table [T} 

bution of the leptons. In the MRSSM, however, the two slepton mass eigenstates are not near one 
another, and so two distinct and well-separated kinematic edges could in principle be extracted. 
This would be a striking signal of slepton flavor violation in the MRSSM. Note also that the electric 
charges of the leptons in this decay are fixed by the conversation of i?-charges. For example, the 
anti-neutralino N!^ can decay into l^^l^, the decay into the same final state for N3 is forbidden. 



5 Discussion 



We have calculated the constraints on /i o e mixing in the MRSSM from the flavor violating 
processes /i — t- 67, — )■ e conversion, and fi — 3e. Given the simplifications stated in Sec. [2| we 
explored LFV in the MRSSM as a function of the parameters Mi, ^(i,m^, and sin2^j- within the 
sub-TeV range. Given the heavier slepton mass set to be m^-^ = 1.5m^-^, we found that the bound 
from — >■ 3e is always less severe than the bounds derived from either /i — )• 67 or — )• e conversion. 
We show the overlapping regions allowed by all constraints in Figs. T5p6 



For right-handed slepton mixing, — t- e conversion in gold nuclei provides the most severe 
constraint - it completely rules out maximal mixing (compare Fig. 15 with Fig. [2]). The situation 



is qualitatively different for left-handed mixing - the most severe bound in this case comes from 
/i — )• 67, as dominant amplitudes (charge radius and magnetic dipole) of /x — )• e conversion interfere 
destructively and opens up a large region in parameter space that satisfies the experimental bounds. 
From Fig. [3] for fj, — t- 67, one sees that maximal mixing is allowed in regions where the bino mass 
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is ~ 500 GeV at = 200 GeV, with a moderate splitting between sleptons. The results suggest 
that the most likely observation of large slepton flavor violation signals at the LHC will occur in 
the left-handed sector. 

Finally, is interesting to consider how the bounds on slepton flavor mixing angles will change 
as the constraints on LFV are strengthened. This is most easily understood by recognizing that 
all of our bounds are proportional to sin^ 20j. In other words, the boundary of the allowed regions 
are contours of constant i?i?bound/ sin^ 29j, where S-Rbound is the bound on the branching ratio of 
a process. In plotting the allowed regions of parameter space in the previous sections of the paper, 
we used of course the current experimental bound. Suppose that in some future experiment the 
bounds are improved, say by a factor of 100. Then, the boundary of the region that satisfy this 
new bound for sin^ 29r = 0.1 is the same as the boundary for the current bound with sin^ 2^,- = 1. 
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Appendix A: Gaugino and Slepton Structure 

To discuss the neutralino masses and interactions more quantitatively, we define the ipB and ipj^^ 

to the be fermion i?-partners of B and H^, respectively. Then we form the Dirac bino and Higgsino 
spinors and their charge conjugates. 



--(S'?)' -'^.KS)' '^-(4)' "^-^(Sf- 

We can also see that the Dirac spinor has an i?-charge —1, whereas A^^ has an i?-charge +1. 
The gaugino mass matrix, Mjv, is shown in the mass term below 

(^-%) ( 'o' ) + (35) 

The mass matrix is diagonalized by a bi-orthogonal transformation; the diagonalized neutralino 
mass matrix, = OIMnOr, obey (M^)2 = OiM^(M;v)^(OL)^ = Or{MnV MNiOnf , where 
0(L,R) are the orthogonal matrices that diagonalize the mass matrix. In this definition, the B and 

content of the i-th neutralino Ni are, Olh and 0^2) respectively. 

We consider mixing between selectrons and smuons only, parameterized as follows: 

h\ ^( cosOj sm9j \ ( ~e\ ^^^^ 



^2 J L,R V -sm^[ cosOi J^^^KfL 



25 



Ml 

^ ^ri (a) = 100 GeV (b) /i^ = 200 GeV 




100 150 200 250 300 350 400 100 150 200 250 300 350 400 

[GeV] [GeV] 



Figure 15: Regions allowed in the parameter space by combining the three constraints for right 
handed sleptons. The constraint from — )• 3e is always less severe than the other two processes in 
the parameter space shown. 
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Figure 16: Same as Fig. 15 but for left handed sleptons. Similar to the right handed case, the 
constraint from — t- 3e is also less severe than the other two processes in the parameter space 
shown. 
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where li represents the sleptons in the mass-eigenstate basis. 

Then slepton flavor violation comes from the interaction terms between a sfermion, fi, a neu- 
tralino, Ni, and a fermion ff. 

- hMuLplOuiGLfip + 0n.2yffRp\) - fRamui^plO luG rJ + ORmffL^]) + h.c. (37) 



where Ul,r are the slepton mixing matrices in Eq. (36). The coupling constants are 



Gl.r = \/29'y,,^„,, and (38) 
V^Mz sm 6u) cos p 

The subscript i on the (s)fermion denotes its generation, subscripts L and R denote the chirality, 
with a and /3 being the flavor indices. The hypercharge of a fermion / is denoted by Yf. From the 
above interaction terms we see that /r and /l have different i?-charges; —1 and +1, respectively. 

The Z-boson only couples to Higgsinos. The ZNN interaction term is 

Z^[7V,7^(Offi20fl,2PL + OLi20Lj2PR)Nj] . (40) 

^ COS t7y; 

One can also write the ZNN coupling in terms of A^'^, 

g 



2 cos 9,1 



Z^,[N-a^{OLi20Lj2PL + ORi20Rj2PR)m], (41) 



Examining the neutralino mixing matrix in Eq. (35), the lightest gaugino receives a negative 
shift, —A < 0, and so the lightest neutralino has mass M^^ = fi^ — A < mci, and thus the lightest 
gaugino is a neutralino. 
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